FTUV/97-80 
IFIC/97-89 
TTP97-26 



OL s {mz) from r decays with matching conditions 

at three loops 



OS 
On 



in 



German Rodrigo^ 1 , Antonio Pich b and Arcadi Santamaria b 

^Inst. fur Theoretische Teilchenphysik, Universitat Karlsruhe 
D-76128 Karlsruhe, Germany 

b Departament de Fisica Tedrica, IFIC, CSIC- Universitat de Valencia, 
E-46100 Burjassot, Valencia, Spain 



> 

o 
o 
t> 

On 

a: 

Oh! 



X 



Abstract 

Using the recent four-loop calculations of the QCD /3-function and the three-loop 
matching coefficients we study the induced errors in a s (mz) obtained from a s (m T ) 
due to the evolution procedure. We show that, when consistent matching and run- 
ning is used at this order, these errors are pushed below 0.0005 in a s (mz)- 



The beta function and the quark mass anomalous dimension govern the evolu- 
tion of the strong coupling constant and the quark masses through the renor- 
malization group (RG) equations, 



da 



(3(a) = -a 2 (/3 + 0ia + 2 a 2 + 3 a 3 ) + 0(a 6 ) 

j, ,, 7m(a) = -a (70 + 7ia + 72« 2 + 73a 3 ) + 0(a 5 ) 
a log u, z K ' 



d log /i 2 
d log fh q 



where a = a s /Tr and rh q is the running mass of the quark q. The coefficients 
of the QCD beta function have been calculated recently in the MS scheme up 
to four loops [1] 
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and also the coefficients of the quark-mass anomalous dimension have been 
calculated at the same order [2,3], 
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Here ( n is the Riemann zeta-function ( ( 2 = k 2 /6, C3 — 1-202056903 . . ., (4 = 
7r 4 /90 and £5 = 1.036927755 . . .) and rif is the number of quark flavours with 
mass lower than the renormalization scale /i. 

Contrary to what happens in momentum-subtraction schemes (MO) this beta 
function and the quark mass anomalous dimension are quark mass indepen- 
dent. Thus, the Appelquist-Carazzone theorem [4], that states that eventu- 
ally the heavy particles decouple at each order in perturbation theory, is not 
realized in a trivial way since coupling constants, beta functions and quark- 
mass anomalous dimensions do not exhibit it. To obtain decoupling in the MS 
scheme we need to build in the decoupling region, ji m, with m the mass of 
the heavy particle, an effective field theory [5] that behaves as if only the light 
degrees of freedom were present. Matching conditions connect the parameters 
of the renormalized low-energy effective Lagrangian with the parameters of the 
full theory. Power suppressed corrections of order 1/m contribute to physical 
observables only through higher order operators but do not affect the match- 
ing conditions for the coupling constant and quark masses. The decoupling of 
the heavy particles is fulfilled in physical quantities once they are expressed 
in terms of the couplings in the effective theory. 

Some time ago it was checked explicitly at three loops [6] that, when the 
appropriate matching conditions are taken into account, the evolution of the 
strong coupling constant from low energies to high energies does not depend on 
the particular choice of the energy scale used to pass a heavy quark threshold. 
The residual dependences that appear in the perturbative calculation are just 
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an estimate of the effects of the higher order corrections. Very recently [7] the 
analysis has been extended to four loops and the appropriate coefficients have 
been computed. 

In this paper we use the recently calculated four-loop MS scheme QCD beta 
function and the quark mass anomalous dimension [1-3] to obtain the logarith- 
mic pieces in the matching conditions from a different point of view. Then, we 
obtain a very convenient analytic form for the running of the QCD coupling 
constant at four loops and compare it with other solutions. Finally, we use 
these results and the non- logarithmic coefficients computed in [7] to analyze 
the impact of the matching conditions on the error induced in a s (mz) if it is 
obtained from a s {m r ) when passing the thresholds of the c and the b quarks. 

Matching conditions in QCD relate the strong coupling constant, a nf , and 
the running mass of the light quarks, fn q , nf) in the full theory with rif flavours 
with the effective strong coupling constant, a n/ _i, and the effective light quark 
masses, m g n/ _ 1 , of the effective theory with rif — 1 flavours through a power 
series in a n/ _i 



O'nf(fJ'th) — a ns ~i{^th 
rh q>nf (fith) =m q ,n f -i(lJ>th) 
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with coefficients that depend on x = log(/^/m 2 ) where m is some RG- 
invariant mass of the heavy quark (for instance the RG-invariant-MS mass, 
rh(rh), or the perturbative pole mass M) that has been integrated out at the 
energy scale fith- In order to obtain a good approximation using only the first 
few terms in the perturbative expansion, we have to evaluate matching condi- 
tions in a region where /it/i/m ~ 0(1). However, the result of these calculations 
should not depend on exactly which \i t h is chosen. 

Note that in contrast to other analysis [7-9] where the effective couplings 
are expressed in terms of the couplings of the full theory we directly write 
the inverted relation since we are interested in the evolution of the QCD 
Lagrangian parameters from low energies to high energies. Note also that in 
order to simplify as much as possible the matching conditions we have taken as 
a reference mass, m, a RG invariant mass instead of the running mass m(n t h) 
evaluated at the threshold scale fj, t h- This makes matching conditions for the 
a s 's independent of the anomalous dimensions. 

The functions C k and H), are, in general, polynomials in x. The coefficients 
multiplying the logarithms of the heavy quark mass are determined just by 
the RG, that is, they are a function of the beta function and the quark mass 
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anomalous dimension of both the effective theory with rif — 1 flavours and the 
full theory with rif flavours. The non- logarithmic coefficients, however, have 
to be be evaluated explicitly for each particular renormalization scheme. 

We apply the renormalization group equations, eq. (f ), to both sides of eq. (4) 
and eq. (5). Identifying order by order in the effective strong coupling constant, 
a nf -i, we obtain for the Ck and the Hk functions a set of coupled first-order 
linear differential equations depending only on the beta and the gamma func- 
tions of the full and the effective theories. By solving them, and using the 
known beta functions we find for the Ck functions 
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where c 2i0 , c 3i0 , c?2,o an d ^3,o are arbitrary constants coming from the integra- 
tion of the differential equations. They depend on the renormalization scheme 
and on the RG-invariant reference mass m chosen. They can be determined 
only by evaluating some Green functions with both the full and the effective 
theories, in a particular mass-independent renormalization scheme, and then 
require they are the same, up to terms 0(l/m), for values of the renormaliza- 
tion scale just around the threshold. Note that, in order to simplify the results, 
we have set c lj0 = d lj0 = 0, which is the MS result with the usual dimensional 
regularization prescription, Tr{I} = 4, with the trace taken in Dirac space. 

If the RG-invariant-MS mass is used as a reference scale, that is m = fh{m) 
the coefficients one obtains are [7,9] 

11 82043 575263 2633 89 

C2 '°-~72' C3 <° " 27648 C3 " 124416 + 31104 n/ ' " _ 432 " (8) 



If the pole mass is used as a reference scale, that is m = M the coefficients 
one obtains are [7] 
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c 2,0 — 
C3,0 = 



(9) 



while ^2,0 does not change, and o^o, when known, has to be shifted by a factor 



Using these coefficients we find complete agreement with ref. [7] also for the 
logarithmic contributions to the matching for both elections of the reference 
scale. Note that the authors of ref. [7] present the inverse relations. 

For consistency, matching conditions at n loops have to be considered together 
with running of the MS parameters at n + 1 loops. The C3 and the if 3 coeffi- 
cients depend on at most the three loop beta function and the three loop quark 
mass anomalous dimension. However, they have to be used together with run- 
ning at four loops, where the recently calculated four loop beta function and 
the quark mass anomalous dimension enter the game. 

Knowing the beta and gamma functions at four loops we could also obtain the 
logarithmic pieces of C 4 and if 4 , however since the coefficients c 4i0 and d^ 
are not known in any particular scheme and since these matching conditions 
should be used together with running at five loops, which is also unknown, we 
do not present them. 

In the following we sketch the formulae used to compute the running of a s 
at four loops, which is the required order if matching conditions are used at 
three loops. 

We solve the full four-loop RG equation for the strong coupling constant as an 
expansion in the solution of the two-loop RG equation. In [6] we obtained an 
expression for the the running of the strong coupling constant as an expansion 
of the strong coupling constant obtained at one loop. We now improve that 
expression by resumming some of the leading dependences proportional to f3\. 
This amounts to expanding around the approximate two-loop solution instead 
of the one-loop solution. 

At the required order we have 



+ 10/27. 



a 



(jm) = a^(jj) (l + c^Xa^/i)) 2 + c 3 (/i)(a( 2 V)) 3 ) , 



(10) 



where aP\n) is the approximate two- loop solution 
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and 



cM = b 2 (l-K) , 

cM = I (l - K 2 ) + hb 2 K (K - 1 - L) + | (L 2 - (1 - ^) 2 ) , (12) 
with K = 1 + p a(/i ) log(/i 2 //ig), L = \ogK and 6 fe = /3 fc /A)- 

Although these expressions are slightly more complicated than the usual ex- 
pansion in 1/ log(AQCT>)[7], they are more convenient because the coupling 
constant at an arbitrary scale is given explicitly in terms of the strong cou- 
pling constant at some reference scale /io, which usually one takes equal to 
mz or to m T . The standard formulae, however, are expressed in terms of 
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Fig. 1. Values one obtains for ah 1 (mz) starting from o^^w-t) = 0.336 at 2,3 and 
4 loops by using the different approximations discussed in the text: i) the usual 
1/ \og(h-QCD) expansion [7] (white), ii) numerical integration of the renormalization 
group equation (soft hatching), iii) our expansion in eq. (10) (hard hatching). 

Aqcd which has to be adjusted to the initial conditions for a s . So, the value 
of A.qcd to be used depends on the number of light flavours considered and 
on the number of loops considered. Moreover it seems that the 1/ log(AQC£>) 
expansion converges more slowly than our expansion. We show in fig. 1 the 
values obtained for af\m z ) starting from af\m T ) = 0.336 by using three 
different solutions of the RG at 2,3 and 4 loops : i) the usual 1/ \og(kqcD) 
expansion [7] (white), ii) numerical integration of the renormalization group 
equation (soft hatching), iii) our expansion in eq. (10) (hard hatching). Clearly 
our expansion gives much closer results to the numerical integration and seems 
to converge faster than the l/log(AQ CD ) expansion. Therefore we will use it 
in the following. 

Let us turn to the phenomenological applications of the matching conditions 
and the four-loop running solutions. 
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The most precise determinations of a s are obtained from hadronic r decays [10] 
at rather low energies and from hadronic Z decays at LEP energies [11] . To 
compare these two results one has to connect the strong coupling constant in a 
theory with three flavours at a scale /x = m T with the strong coupling constant 
in a theory with five flavours at a scale \x — rriz- Therefore, two thresholds have 
to be passed, the threshold of the c-quark and the threshold of the 6-quark. 
Moreover, running between a wide range of scales has to be performed. Note 
that, although m c < m T , results for a s {m T ) are usually presented in a theory 
with only three quark flavours. That is because c-quarks cannot be really 
produced in r decays and they only enter in loops. Therefore it is appropriate 
to use an effective theory in which the c-quark has been integrated out. Power 
corrections of the form m^/m^ can be included in the effective theory and have 
been computed [12,9]. They are very small and they are taken into account in 
the extracted value of a s (m T ). 

Given the present accuracy of both experimental measurements and theoret- 
ical calculations of hadronic decays of the r lepton, it is important to calcu- 
late very precisely the connection between coupling constants when passing 
thresholds and to estimate the remaining errors in the calculation. It is in 
this analysis where the three loop-matching conditions for a s and four-loop 
running studied above are relevant. 

In the following we obtain a s (m z ) by using as starting point a s {m T ) and we 
estimate the residual errors due to the matching conditions and running for 
the different approximations used. We follow the same procedure as in [6]. At 
low energies, /i = m T = 1777.0 ±0.3 MeV we know [10] a s (m T ) = a^(m T ) = 
0.35 ± 0.02. From this we can obtain a^(ii^ h ) at some matching point if th 
around m c by using the renormalization group with rif — 3, then, by using 
eq. (4) with m = rh c and rif = 4 we obtain a^(fx^ h ). Now we use again 
the renormalization group with rif — 4 to obtain aS*\ii b th ) at some matching 
point ix\ h around fhb and use again eq. (4) with m = rhb and rif = 5 to obtain 
°^s\lAh)- Finally we use the renormalization group with rif = 5 to obtain 
af\rriz) = a s (m z ). The final result will depend on the precise values used 
for if th and /i b h and this dependence gives an estimate of the errors which arise 
because the truncation of the perturbative series in the matching conditions. 
In addition, matching conditions also depend on the masses of the quarks, and, 
although they are very well known, their actual value can affect the final result 
for a s {m z )- The induced error due to the uncertainty in the quark masses is 
dominated by the one-loop matching equation. Then, we can estimate this 
error as (Aa s (mz)) / ce s (mz) ~ a s (m q ) / (3n)(Am q ) /m q . We use always as a 
reference scale the RG-invariant-MS mass of the quarks, and therefore the 
coefficients in eq. (8). For the quark masses we take the last values in the 
literature: for the 6-quark mass m b (fh h ) = 4.13 ± 0.06 GeV [13]. For the c- 
quark mass we take rfi c (rh c ) = 1.31 ± 0.06 GeV, [14] (see also [15]). 
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We study the effect of varying the scale at which matching is performed inde- 
pendently for the c and the 6 quarks. 
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Fig. 2. a s (mz) obtained by running the coupling from a s (m T ) = 0.35, as a 
function of the matching point taken to cross the 6-quark threshold. Dotted line: 
two-loop beta functions and one-loop matching conditions. Dashed line: three-loop 
beta functions and two-loop matching conditions. The hatched band is obtained 
with four-loop beta functions and three-loop matching conditions when the 6-quark 
mass is varied within its error interval. 

First we fix [i c th = m c and vary [x\ h in 2 — 20 GeV. Figure 2 shows (dotted line) 
our result for two-loop running (only the first term of eq. (10) is taken into 
account) and one-loop matching conditions (only the C\ coefficient is consid- 
ered). We plot with dashed line the results for three- loop running and two- 
loop matching conditions (two terms in eq. (10) and C2 included in eq. (4)). 
For these two lines we took central values for the 6-quark mass. Finally, the 
hatched area gives the results for four-loop running and three-loop matching 
conditions when the 6-quark mass is varied within its error interval. For the 
central value of the strong coupling constant extracted from tau decays we 
find that varying the 6-quark threshold scale in the range n\ h = 2 — 20 GeV 
two-loop running and one-loop matching conditions induce an error of 0.0006 
on the strong coupling constant at the Z-boson mass scale. With three-loop 
running and two-loop matching conditions the error decreases to 0.0002. For 
four-loop running and three-loop matching conditions we get and error for the 
central value of the 6-quark mass of 0.00009. The uncertainty in the 6-quark 
mass induces in this case an additional error of 0.00003. 

To study the errors induced in passing the c-quark threshold we fix n b th = m b 
and vary \i\ h in the range \i\ h = 1 — AGeV . Then, we find an induced error of 
0.0005, 0.0002 and 0.0001 to each order respectively. The uncertainty in the 
c-quark mass introduces an additional error of 0.0001. 

To analyze the combined effect of passing the two thresholds we have repre- 




8 



sented the value of a s (m z ) obtained at four loops as a function of ii c th and 
H h th as a contour plot in fig. 3 by taking \og(n\ h /\GeV) and \og{nl h /lGeV) 
in the x and y axis respectively. The different contour lines are obtained for 
a s {m z ) = 0.12214 to a s {m z ) = 0.12182 in steps of 0.00004. We see that there 




Log(Mth> 

Fig. 3. a s (mz) as a function of the matching point scales \i1 h and ji\ h . The contours 
are for a s {m z ) = 0.12214 to a s {m z ) = 0.12182 in steps of 0.00004. The area within 
the strait lines is obtained by taking m c < ii c th < n\ h < rn|//i^. 

is a maximum approximately for a s (m z ) = 0.1222. The minimum depends on 
how far away from m& and m c we take the matching scales (/% h and [/% h . We 
choose m c < }/% h < n b th < fhl/[/% h , which are represented by the three strait 
lines in fig. 3. Then we obtain an estimate of the error due to the unknowledge 
of the scales \x\ h and /j^ h of about Aa s (m z ) = 0.0002 by taking half the differ- 
ence between the maximum and the minimum values of a s (m z ). If in addition 
to this we also consider the error induced because the uncertainty in the input 
values of the masses of the c and the b quarks, we obtain that the total error 
in a s (m z ) due to the matching procedure is Aa s (m z ) = 0.0003. Finally we 
can also include the error in the extracted value of a s (m T ), which is the domi- 
nant one. We obtain the values of a s {m z ) starting from a s {m T ) = 0.35 ± 0.02 
and take the average of the maximum and minimum to get the central value 
and half of the difference to obtain the error. This gives Aa s (m z ) = 0.0021. 
Combining lineally all errors we obtain 

a s (m z ) = 0.1219 ±0.0024 (13) 

Clearly, considering matching conditions at three loops and running at four 
loops reduces the errors due to the matching procedure by a factor two and 
renders them below 0.0005 in a s (m z ). This has to be compared with the full 
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error which is of about 0.0025. 



By using the recent four-loop QCD calculations we have obtained the matching 
conditions for a s (/x) and m g (/f) when crossing a quark threshold. Then, we have 
solved the renormalization group equations at four-loops and have obtained a 
very convenient analytic form for the running of the QCD coupling constant. 
We have compared the results obtained by using this solution with the results 
obtained by using the usual l/log(AQcr>) expansion and with the numerical 
solution at the different orders and found that it gives a good approximation 
to the numerical solution and converges faster than the usual l/log(AQC£)) 
expansion. Finally we have used these results to study the effect of the c- and 
the 6-quark thresholds in the evolution of the strong coupling constant from 
fi = m T to fi = rriz- This analysis is very important given the present accuracy 
of both determinations of the QCD coupling constant, a s {m T ) and a s (mz)- 
We found that the total error induced in a s (mz) starting from a s (m T ) due to 
the matching and running procedures is 0.0003 when matching and running 
evolution are used at four loops. 
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